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Band bending inversion in Bi 2 Se3 nanostructures 
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Abstract 

Shubnikov-de-Haas oscillations were studied un¬ 
der high magnetic field in Bi 2 Se 3 nanostructures 
grown by Chemical Vapor Transport, for different 
bulk carrier densities ranging from 3 x 10 19 cm 3 
to 6 x 10 17 cm~ 3 . The contribution of topological 
surface states to electrical transport can be iden¬ 
tified and separated from bulk carriers and mas¬ 
sive two-dimensional electron gas. Band bending 
is investigated, and a crossover from upward to 
downward band bending is found at low bulk den¬ 
sity, as a result of a competition between bulk and 
interface doping. These results highlight the need 
to control electrical doping both in the bulk and at 
interfaces in order to study only topological sur¬ 
face states. 

A well-known limitation to the study of metallic 
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surface states in many 3D topological insulators 
is their finite bulk conductivity 12 . This resid¬ 
ual doping often originates from native defects 
in the crystal. In Bi 2 Se 3 3 4 , electrical doping in¬ 
deed happens through thermodynamically stable 
Se vacancies acting as double donors 56 , induc¬ 
ing a residual bulk carrier density^ 2 5 which can 
be as large as 10 19 cm 3 . Therefore, the full con¬ 
ductance results from the contributions of both 
bulk and topological surface states (TSS). Study¬ 
ing the TSS then requires to identify the their 
own contribution to charge transport. In Bi 2 Se 3 
nanostructures, this was proved possible by study¬ 
ing quantum transport: Aharonov-Bohm effect 7 ®, 
Conductance Fluctuations 9 , and Shubnikov-de- 
Haas oscillations (SdHO) 2 10 19 . So far, the anal¬ 
ysis of SdHO in BHSes remains unclear. Some 
studies reported on single-band SdHOfMEI, 

which 

were attributed to either bulk carriers or TSS. In 
contrast, other studies reported on multi-band 
SdHO 2 16 19 and emphasized the difficulty to 
identify and distinguish between their bulk, TSS, 
or charge-accumulation two-dimensional electron 
gas (2DEG) origin. In this context, an impor¬ 
tant aspect which must be considered is the band 
bending (BB) occuring at the surface/interface 
due to charge transfer 1 . It was investigated by 
ARPES 20 25 and transport measurements 211 ^, 
with seemingly contrasted reports of upward 
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Figure 1: a) Temperature dependence of longitudinal resistance R xx . Insert: SEM and AFM picture of 
sample E, and AFM section along the indicated line, b) Perpendicular field dependence of R xx of sample E; 
inset: second derivative of R xx displaying SdHO. c) Perpendicular field dependence of the Hall resistance 
R xy of sample E («h all = 3.2 ±0.1 x 10 13 cm 2 is the total density determined from high field asymptote 
above 40T). 


(UBB) ™M I or downward (DBB)™!] band 
bending. Yet, this effect is mostly overlooked in 
the analysis of SdHO, and a systematic study of 
both multi-band transport and band bending in 
BESe 3 nanostructures is still missing. 

In this letter, we study Shubnikov-de-Haas os¬ 
cillations in nanostructures of Bi 2 Se 3 under high 
magnetic fields for different bulk carrier densi¬ 
ties ranging over two orders of magnitude, from 
3 x 10 19 cm~ 3 down to 6 x 10 17 cm~ 3 . All con¬ 
tributions to charge transport are separated and 
identified by their temperature and electrical back- 
gate voltage dependence. Band bending is sys¬ 
tematically investigated, and a crossover from 
UBB to DBB is evidenced when bulk doping is 
reduced. Remarkably, a strong DBB induced by 
the SiCB substrate results in the formation of a 
potential well at the interface where a confined 
2DEG coexists with TSS. This study underlines 
the importance of controlling the surfaces quality, 
in addition to solely reducing the bulk residual 
doping, in order to study only the TSS. 

The samples studied in this paper are monocrys¬ 
talline Bi 2 Se 3 nanostructures grown by Chemi¬ 
cal Vapor Transport (for growth details scc 2,s ). 
Cr/Au contacts were prepared using standard e- 
beam lithography and metal lift-off, following an 
in-situ smooth Ar-ion etching to achieve good 
ohmic contacts. All samples show a metallic be¬ 
havior, as expected for degenerate BEScs^ (see 
fig.l.a). Five samples were studied, with different 
bulk doping: high density (A: 1.4 x 10 19 cnV 3 ; B: 


2.6 x 10 19 cm~ 3 ; C: 1.4 x 10 19 cnV 3 ), intermediate 
density (D: 4.3 x 10 18 cm~ 3 ) or low density (E: 

6.6 x 10 17 cm~ 3 ) (see table 1). 



Figure 2: Schematics of BB: left, UBB with A the 
depletion length and AErb > 0 the band bending 
energy; right, DBB. 


For each structure, magneto-transport was stud¬ 
ied down to 4K and up to 55T using a pulsed 
magnetic field, except for structure A which was 
studied under a static field up to 15T. From Hall 
bar measurements, the longitudinal resistance 
Rxx shows a weak amplitude SdHO superim¬ 
posed to a large magnetoresistance background 
(see fig. 1 .b) and a non-linear transverse resis¬ 
tance R xy , which is typical of multiband trans¬ 
port (see fig.l.c). More details on the different 
contributions to electronic transport properties 
can be inferred from SdHO, which quantitative 
analysis is better performed on the second deriva¬ 
tive d 2 R xx /dB 2 . A Fast Fourier Transform (FFT) 
analysis was performed on SdHO to extract the 
frequencies corresponding to the different elec¬ 
tronic populations. In order to identify them, the 
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Figure 3: Magneto-transport study of strongly doped sample A (thickness t=27nm): a) SEM picture of the 
contacted sample A; b) 2 nd derivative of R xx against inverse magnetic field, at two gate voltages; c) FFT 
of SdHO at different gate voltages (curve at V g = -110V shifted for clarity); d) Temperature dependance of 
the amplitude of some SdH peaks (points at 13.2T shifted for clarity); e) Schematic of the band bending 
in sample A. 


temperature and back-gate voltage dependences 
were studied. Indeed, bulk carriers and TSS can 
be told apart through their effective mass m * 2 1 
which is extracted from the temperature depen¬ 
dence of SdHO using the Fifschitz-Kosewitz (FK) 
formula 29 . Moreover, the interface TSS (located 
at the interface between Bi 2 Se 3 and Si 02 ) is gate 
sensitive, contrary to bulk states and surface TSS 
(TSS at the interface with vacuum) 19 . From the 
SdH frequency associated with each population, 
the carrier density, the Fermi energy and the band 
bending AEbb = £p — (F sur f/int — 180meV) are ex¬ 
tracted (see fig.2), assuming that the Dirac point is 
180meV below the bottom of the conduction band 

(CB)SS 

A typical strongly-doped ribbon (sample A) is 
shown in figure 3.a. Above 6 T, the resistance 
exhibits SdHO, that are visible up to 15T (see 
fig.3.b). Those oscillations show a clear multifre¬ 
quency pattern: from ~ 11 T on, the oscillations 


display a third frequency significantly higher than 
below 11T. The FFT reveals three main peaks 
(similarly to earlier results^), at /bi = 165 ± 5T, 
/b 2 = 180 ± 5T and /b 3 = 245 ± 5T, which cor¬ 
respond to three distinct electronic populations 
having different carrier densities. In order to iden¬ 
tify each contribution, SdHO were measured at a 
back-gate voltage of V g = 0 and -110V. The Fourier 
transforms are compared in figure 3.c. As pointed 
out by the dotted lines, both / B 2 and / B 3 do not 
shift at -110V, while /bi is clearly shifted to a 
lower frequency. This indicates that /bi corre¬ 
sponds to the interface TSS. Hence, we calculate: 
Mint = /bi x ( e/h ) = 4.0 ±0,1 x 10 12 cm~ 2 , with 
e the electronic charge and h the Planck constant. 
Knowing the Fermi velocity vp = 5.4 x 10 5 m.s 1 
determined by ARPES 3430 i the interface Fermi 
energy is E- m — 251 ± 3meV and m* nt = E un / vp = 
0.15m e is the associated cyclotron mass. The 
two other frequencies can be separated by their 
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Sample 

w surf 

(cm- 2 ) 

^surf 

(meV) 

AEBB-surf 

(meV) 

«B 

(cm- 3 ) 

E f 

(meV) 

tt int 

(cm- 2 ) 

Ejnt 

(meV) 

AEBB-int 

(meV) 

A 

5.9xl0 12 

305 

+36 UBB 

1.4x10 19 

161 

4xl0 12 

251 

+90 UBB 

E 

2.9xl0 12 

215 

-14DBB 

6.6xl0 lv 

21 

1.3xl0 13 

459 

-258 DBB 

Sample 

^surfl 

Esurl'l 

A^BB-surfl 

riB 

E f 

n surf2 

Esurf2 

AEBB-surf2 

B 

l.OxlO 13 

396 

+35 UBB 

2.6xl0 19 

251 

3.9xl0 12 

247 

+184 UBB 

C 

6.3x10 12 

317 

+38 UBB 

1.5xl0 19 

175 

2.7xl0 12 

206 

+149 UBB 

D 

2.7xl0 12 

207 

+48 UBB 

4.3xl0 18 

75 

1.2xl0 12 

136 

+119 UBB 


Table 1: Table of the carrier densities, band bending and Fermi energies for bulk and TSS, extracted from 
SdHO analysis. Top: devices A and E, for which the interface is identified by its back-gate voltage 
dependence. Bottom: devices B,C,D, measured without an electrical back gate. 


temperature dependences. Since the onset of 
the highest frequency /b 3 is 11T, the tempera¬ 
ture dependence of the SdHO amplitude below 
11T arises only from the two populations asso¬ 
ciated to /bi and /b 2 - Although no quantitative 
value of m* for one single population can be ex¬ 
tracted, LK fits for SdHO peaks below 11T yield 
an averaged effective mass m* — 0.14 ± 0.0 lm e 
(see fig.3.d and Supplementary materials). This 
is consistent with /b 2 originating from bulk carri¬ 
ers. Indeed, = 0.12 — 0.15m e from bulk and 
m* nt from interface are very similar, and tem¬ 
perature dependences are similar for both popu¬ 
lations. Above 11T, m* extracted from LK-fits 
rises to 0.19±0.02m e . This is in agreement with 
the last frequency /b 3 having a higher effective 
mass than the previous two. Therefore we at¬ 
tribute /b 3 to the surface TSS, with an energy 
Esurf = 305 ± 3meV and a corresponding effective 
mass m* urf = 0.185m e , in agreement with the LK 
fits. This is also consistent with the difference 
found in the onsets of the SdHO, which depend 
on the carrier mobility of each band. Indeed, the 
Drude mobility /i = ete/m* depends on the disor¬ 
der strength (characterized by T e = / e /vp, where l e 
is the elastic mean free path) and on the effective 
mass 1 . Here, the carrier density is 50% higher at 
the surface than at the interface. As doping mostly 
originates from double donor Se vacancies 631 , a 
higher doping means a stronger disorder, and a 
reduced Z e . Simultaneously, for TSS, increasing 
energy implies increasing m*. Both effects tend to 
reduce the surface TSS’ mobility compared to the 
bulk and interface one. 

Assuming a quadratic energy dispersion for the 
bulk, the BB can be calculated as shown in fig¬ 


ure 3.e. An UBB is observed in this struc¬ 
ture, as expected for a strongly doped structure 
(«B = 1.4 x 10 19 cnU 3 ). Indeed, if surface dop¬ 
ing only originates from bulk charge transfer and 
not from additional surface doping, the charges 
from the bulk are transferred to the surface over 
a depletion length A, resulting in UBB. Here, 
A > 3 — 4nm assuming a complete depletion of 
the bulk over this range. A can be increased by 
decreasing the bulk carrier density; ultimately, for 
low enough doping and for nanostructures thinner 
than 2A, the bulk could be fully depleted 1 . 

Three other samples (B,C,D), with comparable 
or lower doping, were studied under magnetic 
field up to 55T. As for sample A, the FFTs of the 
multi-frequency SdHO exhibits three peaks, asso¬ 
ciated with bulk and two TSS (see Supplementary 
materials). Because of the absence of a back-gate 
voltage on these structures, surface and interface 
TSS cannot be told apart. Still, all those samples 
show UBB, as can be seen in table 1. For all these 
structures, AErb is of the same order of magni¬ 
tude (~50 to ~150meV); but the effect, small in 
strongly doped structures (for A,B,C: Ef lies in 
the CB at the surfaces), becomes much stronger 
for the intermediate doping (sample D), where 
Ep enters the bulk gap for one surface (E sur f 2 < 
180meV). This confirms the possibility, for thin 
enough structures with doping in the 10 18 cnU 3 
range, to achieve bulk depletion through UBB. 
However, for this doping, the depletion length A 
would be around 3nm for a complete bulk deple¬ 
tion and 6nm for a triangular depletion potential, 
so that bulk depletion can only be achieved for a 
thickness below 12nm, which constitutes a chal- 
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Figure 4: Magneto-transport study of weakly doped sample E (t=79nm): a) SEM picture of the Hall-bar 
designed flake E; b) 2 nd derivative of R xx against inverse magnetic field, at two temperatures to show re¬ 
producibility (curve at -110V in Supplementary materials); c) FFT of the SdHO at different gate voltages; 
d) Temperature dependance of the amplitude of a SdH peak; e) schematic of the band bending in sample 
E. 


lenging task nowadays. To achieve bulk depletion 
over a larger thickness would require to further de¬ 
crease the bulk carrier density. However, this as¬ 
sertion supposes that BB remains upward, which 
is not obvious, as discussed below. 

Magneto-transport results on the weakly-doped 
nanostructure E are presented in figure 4. Com¬ 
plex SdHO are observed in R xx (see figure 4.b), 
with the first oscillations observed at ~6T. The 
multifrequency behavior arises above 15T, which 
imposes to work at high magnetic fields to study 
all populations. The FFT reveals a much more 
complex pattern than for the previous samples 
(see figure 4.c). Five main peaks are identified, at 
/bi = 24 ± 3T, / B2 = 90 ± 4T, / B3 = 121 ± 5T, 
/b 4 = 305 ± 20T and /bs = 550 ± 30T. Impor¬ 
tantly, the lowest frequency fgi corresponds to the 
population with the highest mobility, which gives a 
single-band behavior below 15T. The LK fit of the 


temperature dependence of the SdHO amplitude 
below 15T gives a m* — 0.129±0.01m e (see fig¬ 
ure 4.d), which is consistent with bulk states 32 33 . 
It is also inconsistent with TSS, as for such /bi 
the chemical potential would be 64meV and the 
corresponding cyclotron mass ttij SS — 0.039m e . 
The lowest frequency is therefore unambiguously 
identified as bulk states, with £p = 21meV and 
?7 b = 6.6 x 10 l7 cm 3 . The angle dependence of 
the SdHO, presented in Supplementary materials, 
confirms that all other frequencies than /bi tend 
to disappear when tilting the magnetic field from 
perpendicular to in-plane. 

In order to identify the remaining populations, the 
back-gate voltage dependence of SdHO was stud¬ 
ied up to -32V. The FFT allows us to separate the 
gate-sensitive (/b4,s) and gate-insensitive (/bi,2,3) 
frequencies. Out of the four 2D bands /b2,3,4,5> 
the only consistent model is to associate /b 3 and 
/bs to surface and interface TSS, respectively (see 
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Supplementary materials). The topological carrier 
densities are then n sur f = 2.9 ±0.1 x 10 12 cnU 2 
and n m = 1.33±0.07x 10 13 cm~ 2 , corresponding 
to energies £ sur f = 215 ± 4meV and E- m — 459 ± 
12meV. As a consequence, the strong DBB near 
the interface induces a quantum well, into which a 
confined 2DEG from bulk origin may form. Mod¬ 
eling this well with a triangular potential, one can 
calculate for a thickness A = 9.9 ± 0.9nm that 
two confined sub-bands will be partially filled, 
with energies 40 ± 2 and 143 ± 8meV. They per¬ 
fectly match the remaining /b 2 and /b 4 (40 ± 2 
and 136 ± 8meV) (see fig.4.e). The thickness 
of the quantum well is in agreement with previ¬ 
ously reported DBB 20 23 25 . The total carrier den- 

i =5 

sity, computed as n m — «b x (t — A) ± Y w 9 d = 

i =2 

3.0±0.2 x 10 13 cm 2 taking all five contributions 
into account, is confirmed by Hall measurement 
(^Haii = 3.2±0.1 x 10 13 cm~ 2 , see figurel.b). 
Although the bulk is very weakly doped, BB is 
here downward. This reflects the fact that the in¬ 
terface doping is much higher than the bulk one. 
Such a strong effect on the BB is only visible 
because of the small bulk doping. For a simi¬ 
lar surface doping in sample B (surface 1) but 
with a much larger bulk carrier density, BB re¬ 
mains upward due to the large Fermi energy in 
the bulk. This points at the necessity to control 
both the residual bulk doping and the additional 
surface/interface doping in order to produce bulk- 
compensated structures. In particular, the very 
strong DBB near the interface suggests that a Si02 
substrate favors a strong interface doping, which 
could result either from a locally-increased disor¬ 
der in Bi 2 Se 3 or from chemical bonds. The use of 
an alternative oxygen-free substrate (GaAs, SiC, 
GaN...) could maybe overcome this effect. 

The difference in the observed mobilities is also 
consistent with the type of band bending and the 
assumption that disorder is stronger at interfaces 
with respect to the bulk. In sample E, the situation 
Mint < M surf < Mb is well understood by Z e . b > / e . T ss 
and m* nt > m* urf ~ m b , as a result of DBB. Indeed, 
bulk SdHO start first (~ 6 T), whereas surface and 
interface SdHO begin at higher fields (~ 147’ and 
~ 27 T, respectively, in the ratio of their cyclotron 
masses). On the contrary, for samples A,B,C,D, 
bulk states and TSS have a more similar mobility, 


which is in agreement with a reduced cyclotron 
mass in the case of UBB. 

In conclusion, we investigated charge transport 
in Bi 2 Se 3 nanostructures grown by CVT through 
a careful analysis of SdHO under high magnetic 
fields. This gives the possibility to identify all 
electronic populations, as bulk states and TSS, or 
even as a massive 2DEG for low bulk doping. We 
confirm that UBB is the usual case for strongly 
doped structures, due to bulk-dominated charge 
transfer at interfaces. Remarkably, DBB is found 
for a rather small bulk doping, which reveals the 
importance of interfaces doping and the influence 
of the substrate. Although UBB could indeed lead 
to bulk-depleted structures, our results emphasize 
the importance to control the interfaces quality so 
as to avoid downward band bending, which is nec¬ 
essary for further investigation and engineering of 
TSS transport. A possible direction maybe could 
be the growth on oxygen-free substrates. 
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